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Abstract 

The problem of fermions in 1 + 1 dimensions in the presence of a pseudoscalar Coulomb potential plus a mixing of 
vector and scalar Coulomb potentials which have equal or opposite signs is investigated. We explore all the possible 
signs of the potentials and discuss their bound-state solutions for fermions and antifermions. We show the relation 
between spin and pseudospin symmetries by means of charge-conjugation and y 5 chiral transformations. The cases of 
pure pseudoscalar and mixed vector-scalar potentials, already analyzed in previous works, are obtained as particular 
cases. The results presented can be extended to 3+1 dimensions. 
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1. Introduction 

The potential generated by a point charge, the Coulomb potential, depends on the dimensionality of space-time 
(see, e.g. [1]). The (l + l)-dimensional Coulomb potential is linear and so it provides a constant electric field always 
pointing to, or from, the point charge. This problem is related to the confinement of fermions in the Schwinger 
and in the massive Schwinger models [2, 3] as well as in the Thirring-Schwinger model [4]. Considered as the time 
component of a Lorentz vector and due to the tunneling effect (Klein's paradox), there are no bound states for this 
kind of potential regardless of the strength of the potential [5, 6]. The linear potential, considered as a Lorentz scalar, 
is also related to the quarkonium model in 1 + 1 dimensions [7, 8]. Although it was incorrectly concluded that even 
in this case there is just one bound state [9], later on the proper solutions for this last problem were found [10-12]. 
However, it is well known from the quarkonium phenomenology in the real 3+1 dimensional world that the best fit for 
meson spectroscopy is found for a convenient mixture of vector and scalar potentials put by hand in the equations (see, 
e.g., [13]). The same can be said about the treatment of the nuclear phenomena describing the influence of the nuclear 
medium on the nucleons [14]. The mixed vector-scalar potential has also been analyzed with the Dirac equation in 
1 + 1 dimensions for a linear potential [15] as well as for a general potential which goes to infinity as \x\ —> oo [16]. 
In both of those last references it has been concluded that there is confinement if the scalar coupling is of sufficient 
intensity compared to the vector coupling. The problem has also been analyzed for pseudoscalar couplings [17]. 

The case in which the mean field is composed of a vector (V,) and a scalar (VA potentials, with V, = -V s , is 
particularly relevant in nuclear physics, because it is usually pointed out as a necessary condition for occurrence of 
pseudospin symmetry in nuclei [18-28]. Furthermore, with an appropriate choice of signs, potentials fulfilling the 
relations V s = ±V, are able to bind either fermions or antifermions [29-33]. 
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Closely related to this is the fact that, in the nucleus, the charge-conjugation transformation relates the spin sym¬ 
metry of the negative bound-state solutions (antinucleons) to the pseudospin symmetry of the positive bound-state 
solutions (nucleons) [30, 34]. Therefore, we believe that this connection between spin and pseudospin symmetry 
obtained by charge conjugation deserves to be more explored. 

The main motivation of this article is illustrate the relation between spin and pseudospin symmetries using charge- 
conjugation and chiral transformations in the Dirac equation in 1 + 1 dimensions, which was already uncovered in 
ref. [35], for the problem of mixed scalar-vector-pseudoscalar Coulomb potential in 1 + 1 dimensions. We thus take 
advantage of the simplicity of the lowest dimensionality of the space-time as much as was done for the harmonic 
oscillator potential [35] and Poschl-Teller potential [36]. This approach is equivalent to consider fermions in 3+1 
dimensions that are restricted to move in one direction [37]. We explore the spectra when it is possible to obtain 
analytical solutions, i.e., in the particular cases when A = V, - V s — 0 or Z = V t + V s — 0. We explore all the possible 
signs of the potentials, thus paying attention to bound states of fermions and antifermions as well. We compare 
both cases A = 0 and Z = 0 to establish the charge-conjugation connection discussed above in the presence of the 
pseudoscalar term. We also discuss the connection between spin and pseudospin symmetries by means of the chiral 
transformation. 

2. The Dirac equation in 1+1 dimensions 

The 1 + 1 dimensional time-independent Dirac equation for a fermion of rest mass m under the action of vector (V,), 
scalar (fA) and pseudoscalar (V p ) potentials can be written, in terms of the combinations Z = V, + V s and A = V t - V s , 
as 

Hi//= Eifi, H = ccT\p + crpnc 2 H-——-Z H-^-A + tT 2 V p , (1) 

where E is the energy of the fermion, c is the velocity of light and p is the momentum operator. The matrices cn, 
cr 2 and cr 3 denote the Pauli matrices, and I denotes the 2x2 unit matrix. The positive definite function <//| 2 = if f, 
satisfying a continuity equation, is interpreted as a position probability density and its norm is a constant of motion. 
This interpretation is completely satisfactory for single-particle states [38]. 

2.1. Parity, charge conjugation and chiral transformation 

The Dirac equation is covariant under jc —> —jc if V p changes sign whereas V s and V, remain the same. This is 
because the parity operator P = exp {is) Pocr 2 , where s is a constant phase and Pq changes x into —x, changes the sign 
of cr ! and 0-2 but not of <x 3 . 

The charge-conjugation operation is accomplished by the transformation if/ c = oqi/A and the Dirac equation be¬ 
comes H c if/ c = -Ei// C , with 

1 I + 0"3 I — (7"3 

H c = ccrip + cr 3 mc 2 -— A-— z + cr 2 V p . (2) 

One see that the charge-conjugation operation changes the sign of the energy and of the potentials V, and V p . In turn, 
this means that Z turns into -A and A into -Z. Therefore, to be invariant under charge conjugation, the Hamiltonian 
must contain only a scalar potential. 

The chiral operator for a Dirac spinor is the matrix y 5 = cr\. Under the discrete chiral transformation the spinor 
is transformed as - y 5 i// and the transformed Hamiltonian H x — y 5 Hy 5 reads 

tt 2 I + °"3 , 1 ~ o"3 v 

H x = ccrip - cr 3 mc + —-—A+—-—Z + o- 2 V p . (3) 

This means that the chiral transformation changes the sign of the mass and of the scalar and pseudoscalar potentials, 
thus turning Z into A and vice versa. A chiral invariant Hamiltonian needs to have zero mass and V s and V p zero 
everywhere. 
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2.2. Equations of motion and the Stunn-Liouville problem 

If we now write the spinor if in terms of its components if T = (if+, iff), the Dirac equation gives rise to two 
coupled first-order equations for the upper, if + and the lower, if components of the spinor: 

—ihcif'_ + mc 2 if+ + Xt f + — iV p if- = Eif + (4) 


-ihci/j' + - mc 2 if- + Aif- + iV p if+ = Eif- , (5) 

where the prime denotes differentiation with respect to x. In terms of if , and if- the spinor is normalized as 
XT dx(\if + \ 2 + \if-\ 2 ) = 1, so that <// + . and if are square integrable functions. It is clear from the pair of coupled 
first-order differential equations (4) and (5) that if+ and if have opposite parities if the Dirac equation is covariant 
under x —> -x. 

For A = 0 with E 4 -me 2 , the Dirac equation becomes 


ftcif' + — Vpifr 

= -i - 7 -- 2~ 

E + me- 


+ 


( 6 ) 


n 2 (E + me 2 )!, + V 2 n + hcV' 

_ _ _ ,// ____ P _P 

2m + 2 me 2 

and for 1,-0 with E 4 me 2 , the Dirac equation becomes 




E 2 - m 2 c 4 
2 me 2 ^ + ’ 




hcf_ + Vpif- 
E — me 2 


(7) 

( 8 ) 


ft 2 (E - me 2 )A + V 2 - hcV’ 

-7— f" + - 

2m 2mc z 


E 2 - m 2 c 4 
2 mc 2 ^ 


(9) 


Either for A = 0 with £ 4 -me 2 or X = 0 with E 4 me 2 the solution of the relativistic problem is mapped into a 
Sturm-Liouville problem in such a way that the solution can be found by solving a Schrodinger-like problem. In the 
case of invariance under reflection through the origin (x —> -x), wave functions with well-defined parities can be built. 
Thus, it suffices to study only the positive half-line and impose adequate boundary conditions on if ., or if- at the origin 
and at infinity. Normalizability demands that if ± —> 0 as x —> oo. Eigenfunctions and their first derivatives which are 
continuous on the whole line with well-defined parities can be constructed by taking symmetric and antisymmetric 
linear combinations of if ± defined on the positive side of the x-axis. 

The solutions for A = 0 with E = -me 2 and I = 0 with E = me 2 , excluded from the Sturm-Liouville problem, 
can be obtained directly from the original first-order equations (4) and (5). They are 


* + =*®e*p[ + J>^] 

( 10 ) 

+ + 2mc 2 )if + 

for A = 0 with E — -me 2 , and 

*_=^exp[-f*^] 

( 11 ) 

<A'+ - = ~/t( A ~ 2 m c 2 )if- 

for E = 0 with E = me 2 , where i/p 1 and if m> are normalization constants. 
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3. The Coulomb potential 

Let us consider 


Y. — k\\x\, A = 0, Vp — kix. (12) 

As we have seen above, the space inversion does not change Z and A but changes V p by —V p . The chiral trans¬ 
formation performs the changes A —> Z, Z —> A, m — > —m, and V p —> —V p . Moreover, because y 5 interchanges the 
upper and lower components, the resulting pair of transformed equations of motion are formally the same, so that 
their solutions have the same energy eigenvalues. This symmetry can be clearly seen from the two equation pairs 
(6)-(7) and (8)-(9), as well as from the isolated solutions (10) and (11), which are converted into each other by this 
kind of transformation. This transformation provides a simple mechanism by which one can go from the results from 
the Z = k\\x\, A = 0, V p = k 2 X case to the case when A = k\\x\, Z = 0, V p — by just changing the sign of m in the 
combinations (E ± me 2 ) and of kn in the relevant expressions. Notice, however, that in this latter case the potentials 
are different and so the eigenenergies will also be different. 

The Dirac spinor corresponding to the isolated solution with E = -me 2 is obtained from Eq. (10). Only for E > 0 
there is a normalizable Dirac spinor, the upper component vanishes, whereas the lower component is an even-parity 
function given by iff- = if/°_e\p [-k 2 X 2 /( 2 fic)J. For E + -me 2 , Eq. (7) takes the form 


2 m 


if/" + I —A x 2 + B\x\ I 


E 2 


2 4 

■ me 


■ tic ki 


2 me 2 


i/s + . 


where 


A = 


r*2 ’ 


B = 


k\ 


-(E + me). 


(13) 


(14) 


As referred above, when A = k\ |x|, Z = 0, V p - kix, one obtains the equation for if/ easily from (9), so that 
one has 


7r 


2m 


- — if/" + -A' x 2 + B'\x\ iff- = 


nrc 4 + Tic k2 
2 me 2 


•If-, 


where 


k 2 

A' = 2 


. 2 ’ 


B' = 


k\ 


■(E - me 2 ). 


(15) 


(16) 


The equations presented are similar to the Dirac equations in 3+1 dimensions, with spin and pseudospin symme¬ 
tries, either with one-dimensional or radial linear confining potentials, in this last case for s (l = 0 ) states. 


3.1. Pure scalar and vector couplings 

In this section we present the solutions of eq. (13) when &2 = 0 (A = 0). As already mentioned, the solutions 
of this case and for Z = 0 were already presented in [15]. Here we review those solutions in a different perspective, 
pointing out the relationship between the A = 0 and Z = 0 solutions and also discussing the non-relativistic limits. 
This sets a convenient framework for the discussion of the new results presented in the next subsection. 

When k 2 = 0 (A = 0), (13) reduces to 


Tr 

-—iff” + B\x\if/+ 
lm 


E 2 - m 2 c 4 
2 me 2 


(17) 


For this class of potentials, the existence of bound-state solutions requires B > 0 and thus, from (14), one must have 
E % -me 2 for k\ % 0. Eq. (17) for x > 0 turns into the Airy differential equation 


d 2 if/+ 

dz 2 


- zif/+ = 0 , 


( 18 ) 


where z = ax + b, a - [k\l(Trc 2 ) {E + me 2 )] 1 / 3 , b = —a{E — mc 2 )/k\, which has square-integrable solutions expressed 
in terms of the Airy functions [39]: if/+(z) = CAi(z), where C is a normalization constant. Continuity of i f/+ and 
its derivative at x = 0 imply that the homogeneous Dirichlet boundary condition (if/( 0) = 0) must be satisfied for 
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odd-parity solutions whereas the homogeneous Neumann condition ( di///dx\ x= Q = 0) must be satisfied for even-parity 
solutions, i.e., Ai (b) — 0 for odd-parity and Ai '(b) = 0 for even-parity solutions. 

These quantization conditions have solutions only for b < 0 and can be found numerically. From the definition of 
b one can see that b < 0 always corresponds to \E\ > me 2 regardless of the sign of k\. From the roots of Ai (b) and 
Ai'(fi) we obtain the possible energies as the solutions of a fourth-degree algebraic equation: 

(E - mc 2 )\E + me 2 ) - (heki) 2 \b\ 3 = E 4 - 2mc 2 E 3 + 2m 3 c 6 E - [m 4 c s + {hck x ) 2 |fi| 3 ] = 0 (19) 

which can be recast in the form (with E — me 2 + mc 2 d +) 

Sl(6 + +2)-D = 6 4 + + 2di-D = 0 (20) 

where D = [Hk\/{m 2 c 3 )] 2 {-b ) 3 = K 2 {-b ) 3 and k\ - k\ A c /{me 2 ) 2 = kih/(m 2 c 3 (, A c = H/(mc) being the reduced 
Compton wavelength and k\ the dimensionless strength of the potential. Since D is positive, there is only one positive 
solution according to Descartes' rule of signs for the roots of polynomials (see, e.g., [40, 41]) and so this is the only 
solution for k\ > 0. For k\ < 0, since one have E < me 2 , we can set E = -me 2 - me 2 5 _ and find 

5-{8- + 2 ) 3 - D = 8t + 68l + 128 2 + 8d_ - D - 0 . (21) 

Again, there is only a positive root, and therefore the only root. It is interesting to note that this result is true whatever 
the values of the fermion masses and the coupling constant. Fig. 1 shows the behavior of the energies for the four 



Figure 1: First four energy levels as a functions of the dimensionless coupling constant k\ when ki = 0 and A = 0. The black circles represent the 
non-bound solutions E = ±mc . 


lowest energy levels as a function of k\, when A = 0. From Fig. 1 one sees that all the energies levels for k\ (k\) > 0 
emerge from me 2 and one finds that the lowest quantum numbers correspond to the lowest energies, as it should be for 
particle energy levels. For k\ (k\) <0 the spectrum presents a similar behavior but all the energies levels emerge from 
- me 2 and the highest energy levels are labeled by the lowest quantum numbers. These energy levels can be identified 
with antiparticle levels. This conclusion confirms what has already been analyzed in [42, 43]: the spectrum contains 
either particle-energy levels or antiparticle-energy levels depending on the sign of the coupling constant, but not both 
kind of levels for a certain value of coupling constant, as is the case in the 3+1 Coulomb problem [31]. The same 
behavior occurs for harmonic oscillator potentials as reported in [35]. 

If one takes the non-relativistic limit by setting |/q | <$; 1, then, from (20), one gets for positive k\ ( 5 + <s 1) 


2(5 3 = D 



( 22 ) 


5 





and for negative k\ (<5_ «c 1), from(21), 


( 23 ) 


D 

8 8- = D => 6- = — . 

8 

In terms of k\, the result for 6 (negative energy, anti-fermions) is of higher order (3 times as much) than 6+ (positive 
energy, fermions). From this we may state that there is no non-relativistic limit for negative energy solutions, since 
the solution corresponds to a higher-order term in a l/(me 2 ) expansion. Furthermore, one can check that the solution 
of the one-dimensional Schrodinger equation with a linear confining potential depends on the strength of the potential 
raised to a power of 2/3, as is the case for 6+ in (22). 

In the case when X = 0, i.e., setting A' = 0 and ki — 0 in eq. (15), one gets 


n 2 

+ B ' |j#_ 


E 2 - m 2 c 4 
2 me 2 ^ 


d 2 >//- 

dz' 2 


- z'i/j- = 0 when x > 0 , 


(24) 


where z! = a!x + b ', a' = [k\/(ti 2 c 2 ) ( E - mc 2 )] 1 ^ 3 , b' = -a' (E + mc 2 )/k\. One finds again that one must have E > me 2 
or E < -me 2 when k\ > 0 or k\ < 0, respectively. Also one may notice immediately that a' and b' could be obtained 
from a and b by changing m —> —m, as referred before. So it is not surprising that the solutions of this equation are 
again Airy functions, the values of b' being its zeros or the zeros of its derivative, depending on to the parity of <// . 
Therefore, one can set b' - b. The eigenvalue equation, with k\ > 0, E - me 2 + mc 2 d' + is 

d;((5; + 2) 3 - D = 5'\ + 6 S’l + 12 6'1 + 8 6' + -D = 0 , (25) 


where D is given as before. Eq. (25) is identical to (21) so one can state that there exists one solution with E > me 2 
and that there is no non-relativistic limit for this solution. For k\ < 0, i.e., for E - -me 2 - mc 2 8'_ < 0, one gets an 
equation for <51 identical to (20) for <5+. 



Figure 2: First four energy levels as a functions of k\ when hi = 0 and S = 0. The black circles represent the non-bound solutions E = ±mc 2 


This can be seen from Fig. 2 which shows the behavior of the energies for the four lowest energy levels as a 
function of /q, when X = 0. Thus we can conclude that the solutions for positive and negative energy when X = 0 are 
reversed with respect to the solutions when A = 0 in the sense that, for the same value of |/q[, the energies \E\ are the 
same for positive (negative) k\ for A = 0 as for negative (positive) k\ for X = 0. Probably the most important result 
is that one may regard the positive solutions for X = 0 as intrinsically relativistic, since there is no non-relativistic 
limit. This was already pointed out in ref. [35] quite generally for any potential A when X = 0. One has also shown 
that, as in [35] for harmonic oscillator potentials, one is able to find bound solutions when X = 0 for linear confining 
potentials as is the present case. This agrees with the general finding that in the 3+1 spherically symmetric case there 
are only bound solutions for X = 0 when the vector and scalar potentials are confining potentials [44]. 
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3.2. The addition of a pseudoscalar coupling 

For k 2 4 0 (A 4 0), the equation (13) (A = 0) can be cast into the form 


where 






1 / E 2 - m 2 c 4 — h ck-2 B 2 \ 1 fm 

2 y 2 me 2 2A j h \ A 

E 2 - m 2 c 4 1 tk 2 \ k] (E + me 2 ) 2 

2hc\k 2 \ ~ 2 [ M + 1 ) + jfep 8 Tc 


( 26 ) 


(27) 


The existence of bound solutions of equation (26) is guaranteed by the fact that A > 0, independently of the sign 
of k 2 . Let us define for positive x 


y = ax + /%, 



(h 


2\k 2 \ k\(E + me 2 ) 


he 


2 k\ 


so that (26) transmutes into 


dy 2 


+ 




= 0 . 


(28) 


(29) 


The boundary condition i//+ —> 0 as x —> oo implies that we must seek solutions of (29) which vanish as y —> oo. The 
particular solution of (29) which vanishes for very large positive values of y is called a parabolic cylinder function, it 
is denoted by D v (y). Now, the boundary conditions at x = 0 lead to the quantization conditions 


D v (J3f) = 0 for odd-parity solutions 
Dy (Pf) = 0 for even-parity solutions. 


(30) 


These equations are non-explicit equations for the energy E, since the dependence on E comes about not only through 
Py. but also through v. From the properties of the real zeros of the the parabolic cylinder functions [45] one must have 
v > 0 for D v to have real zeros, and their derivatives D' v have real zeros for v > -0.21, as checked with the program 
MATHEMATICA. If v = 0 D' v has a zero at the origin. This means that the following condition for the energies and 
strengths of the potentials must be satisfied 


E 2 -m 2 c 4 1 / k 2 \ k\ (E + mc 2 ) 2 

2hc\k 2 \ ~ 2 + 1 ) + fep She 


(31) 


When v > 0, solutions with both parities exist. Also, depending on the particular value of v, there are both positive 
and negative zeros of D v and D' v , so that, from (28), there can exist both positive and negative energy solutions for a 
particular pair of {k\,k 2 ) values. As mentioned before, in a purely confining scalar and/or vector potential this is not 
possible. 

The fact that the zeros of D v [45] and also D' v are bounded by the value of v further restricts the existence of 
solutions in this case. This is because the values of /L and v, for fixed k\ and k 2 , are themselves connected by their 
expressions in eqs. (28) and (31) by elimination of the energy E. Numerical tests show that the ratio k\/k 2 cannot be 
big in order to have solutions for this problem. 

The special case of k\ — 0 is worth mentioning. In this case, Py - 0 and v = n - 0, 1, 2, ... is an integer. We have 
then D n (y) = 2~ n/1 e~i yl H n (y/ V2), where H n are the Hermite polynomials of degree n [45]. This corresponds to the 
pseudoscalar Coulomb potential, already analyzed in [17] or, viewed in another way, the so-called Dirac oscillator in 
1 + 1 dimensions studied in [35]. The spectrum has both positive and negative energies for a particular value of k 2 . 
Defining k 2 = k 2 /(tim 2 c i ) as the dimensionless strength of the pseudoscalar coupling, one has explicitly 

E 2 = m 2 c 4 [ 1 + k 2 + \K 2 \(2n + 1)] . (32) 
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We used this solution as a starting point for finding solutions for this problem, fixing K2 and allowing k\ to vary. 
We chose a value of K2 large enough in order to keep the ratio k\ /kj sufficiently small for the range of k\ values 
considered, since, as referred before, otherwise one cannot get solutions to this problem. Because of this restriction, 
the complementary study (fixing k\ and allowing ko to vary, starting from the K 2 - 0 solution) is not feasible. Note 
that the quantum number n in (32) defines the parity of the solutions of the Dirac oscillator and, by continuity, the 
parity of the solutions in the present case. For even or odd n one has even or odd solutions, respectively. Furthermore, 
for K2 > 0, n takes the values 0,1,..but for K2 < 0 the value n - 0 is excluded [35]. In Figs. 3 and 4 we plot 
the solutions found numerically for kt = 5 and Dirac oscillator quantum number n = 0,1,2,3 and for K 2 = -5 and 
n = 1,2,3, respectively. 



Figure 3: First four energy levels as a functions of k\ when ki = 5 and A = 0. 



Figure 4: First three energy levels as a functions of k\ when K2 = -5 and A = 0. 


Now, let us present, as an example of the rule stated at the beginning of the section, the quantization conditions 
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for A = k\\x\, £ = 0 ,E± me 2 and V p = k 2 x. These are 


D p = 0 for odd-parity solutions 

D' fl (fj\) = 0 for even-parity solutions, 


(33) 


where 


and /.(is given by 


Pa 


2\k 2 \ k\(E — me 2 ) 
~hc 2 k 2 


E 1 - m 2 c 4 ! 1 / k 2 \ k\ (E- me 2 ) 2 
2tic\k 2 \ + 2\M~ l ) + \k£ &tic 


(34) 


(35) 


The second-order equation to solve in this case (£ = 0) is the one for the lower component, Eq. (9), because the 
chiral transformation interchanges the upper and lower components. The upper component can be obtained from the 
corresponding first-order equation Eq. (8). Note that for massless particles, the cases A = 0 and I = 0 have the same 
spectrum with the sign of k 2 reversed. 



Figure 5: First four energy levels as a functions of k\ when K2 = -5 and E = 0. 


In Figure 5 we show the solutions for £ = 0 found numerically for k 2 = -5 and Dirac oscillator quantum number 
n — 0,1,2,3 (note that, when solving for £ = 0 and taking the limit k\ —> 0 in the expression for the energy of Dirac 
oscillator (32), the sign of k 2 is reversed). 

Comparing with Figure 3, we see that the plots are identical if we reverse both the vertical and horizontal axes, i.e., 
we get the same solutions as in the case of k 2 = 5 and A = 0 if we reverse the sign of the energy and of k\ , respectively. 
This is because /A. turns into fi\ if we make the changes k\ —> —k\ and E —» —E and, on the other hand, v turns into 
fi when k 2 —> -k 2 and E —> —E, while it is left unchanged by the change of the sign of k\. This illustrates the fact 
that the case £ = 0 can also be obtained from the A = 0 case by applying the charge-conjugation transformation. We 
recall that this transformation performs the changes E —> —E, A —> -£ V p — > -V p , i.e., the changes k\ —> —k\ and 
k 2 —> — k 2 . 

Thus we see that, by adding the pseudoscalar coupling, one allowed for the possibility of having bounded states 
with both positive- and negative-energy solutions in a system with either A = 0 or £ = 0, which can be relevant in 
strong interaction physics where these linear confining potentials can appear. 
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4. Conclusions 


In this paper we have computed and described in detail the bound-state solutions of the 1 + 1 Dirac equation with 
a Coulomb potential with the most general Lorentz structure in spin and pseudospin symmetry conditions. We found 
that there exist bounded solutions for both particles and antiparticles but not for all the range of the coupling strengths. 

At the same time, this work illustrates some general conclusions drawn in previous works about spin and pseu¬ 
dospin symmetry, namely that one can obtain the solutions for 2 = 0 from the A = 0 case, using the charge-conjugation 
and chiral transformations. 

Another important conclusion is that, as previously shown for harmonic oscillator potentials, one is able to get 
bound states in pseudospin symmetry conditions (2 = 0, with or without a pseudoscalar potential), which seems 
to confirm the assertion, made in the context of the 3+1 Dirac equation with radial potentials, that the asymptotic 
behavior of the potentials (confining in this case) is a crucial feature that allows to have such bound states. 

To finish with, one should also refer that these findings may be applied to the one-dimensional linear potential 
problem in a 3+1 Dirac equation and in principle in the spherically symmetric problem also in 3+1 dimensions, 
provided one has £ - 0 or £ = 0 in the spin symmetric (A = 0) or the pseudospin symmetric (2 = 0) cases, 
respectively, £ and £ being the orbital and pseudospin orbital angular momentum quantum numbers, respectively [46]. 
As such, these results might be of relevance to quarkonium phenomenology in 3+1 dimensions. 
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